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Abstract. — In this paper we prove an abstract KAM theorem for infi- 
nite dimensional Hamiltonians systems. This result extends previous works 
of S.B. Kuksin and J. Poschel and uses recent techniques of H. Eliasson and 
S.B. Kuksin. As an application we show that some ID nonlinear Schrodinger 
equations with harmonic potential admits many quasi-periodic solutions. In a 
second application we prove the reducibility of the ID Schrodinger equations 
with the harmonic potential and a quasi periodic in time potential. 
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1. Introduction 

Let ^ : N — > [0, +oo[ so that ^(j) > j for all j > 1. We consider the 
(complex) Hilbert space £|, defined by the norm 

We define the symplectic phase space V = as 
(1.1) V = T' X i% ^ 1%, 

equipped with the canonic symplectic structure: 

n 

dOj A dyj + duj A dvj . 
3=1 j>i 
For (0, y, n, f ) G "P we introduce the following Hamiltonian in normal form 

n 

(1-2) N = Y, ^jiOVj + 2 E + ^D' 

where ^ G M" is an external parameter. 

In [lOj . (see also [llj and a slightly generalised version in [16]) S.B. Kuksin has 
shown the persistence of n— dimensional tori for the perturbed Hamiltonians 
H = N + P with general conditions on the frequencies loj , ilj and perturbation 
P which essentially are the following : Firstly the frequencies satisfy some 
Melnikov conditions and the external frequencies have to be well separated 
in the sense that there exists d > 1 so that roughly speaking (see Assumption 
El below) 

(1.3) %(e)«/. 

Denote by V^'P the phase space given by the weight ^{j) = jP/'^e"'^ where p > 
and a > 0. Secondly, the perturbation is real analytic and the corresponding 
Hamiltonian vector field is so that 

{p > p for d > 1, 
p > p for d = 1, 

where d is the constant which appears in (|1.3|) . For instance, the Schrodinger 
and the wave equation on [0, vr] with Dirichlet boundary conditions satisfy the 
previous conditions, see respectively the KAM results of Kuksin-Poschel [13] 
and Poschel |18] . Indeed the result in |13j is stronger because there is no 
external parameter ^ in the equation. 

Now, if we consider the nonlinear harmonic oscillator 

(1.5) idtu = -d^u + x^u + V{x)u + \u\^'^u, (t,x)GMxM, 
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with real and bounded potential V, we have 0,j ~ 2j + 1, hence d = 1 but the 
Hamiltonian perturbation which is here 

(1.6) P = [ (nn)"+Mx, 

Jr 

does not satisfy the strict smoothing condition (ll.4p (see Section [6] for more 
details). 

The aim of this paper is to prove a KAM theorem (Theorem 12. 3p in the 
case d = 1 and p = p in ()1.3p and ()1.4p . To compensate the lack of smoothing 
effect of Xp we need some additional conditions (see Assumption S]) on the 
decay of the P derivatives (in the spirit of the so-called Toplitz-Lipschitz 
condition used by Eliasson &; Kuksin in [6j) which will be satisfied by the 
perturbation ()1.6p . The general strategy is explained with more details in 
Section 12.31 

Notice that S.B. Kuksin has already considered in [llj the harmonic oscillator 

with a smoothing nonlinearity of type P = ip{\u -k ^\)dx where ^ is a fixed 

Jm. 

smooth function. 



We present two applications of our abstract result concerning the harmonic 
oscillator T = —d'^ + x^. Let p>2 and denote by l"^ the space l\ with ^'(j) = 
jP/2_ 'pj^g operator T has eigenfunctions {hj)j>i (the Hermite functions) which 
satisfy Thj = (2j — l)/ij, j > 1 and form a Hilbertian basis of L^(M). Let 
u = X]j>i '^jhj ^ typical element of L^(M). Then {uj)j>i G ip if and only 
iiu^W := D{TpI'^) = {n G ^^(R) | TpI'^u € L2(R)}. Indeed W is a Sobolev 
space based on T and we can check that 

UP = D{TP''^) = {n G L2(R) I x'^d^u G ^^(R) for a + /3 < p}. 

In this context, we are able to apply our KAM result to ()1.5p and we obtain 
(see Theorem 16.61 for a more precise statement) 

Theorem 1.1. — Let m > 1 be an integer. For typical potential V and for 
e > small enough, the nonlinear Schrddinger equation 

(1.7) idtu = -dlu + x^u + £V{x)u ± e\u\^'^u 

has many quasi-periodic solutions in 71°° . 

Here the notion of "typical potential" is vague. This means that there 
exists rather a large class of perturbations of the harmonic oscillator so that 
the result of Theorem 11.11 holds true (unfortunately our result does not cover 
the case V = 0). Since the definition of this class is technical, we postpone it 
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to Section [6l 

The physical motivation for considering equation ()1.7p (for V = 0) comes 
from the Gross-Pitaevski equation used in the study of Bose-Einstein conden- 
sation (see |15j ). The harmonic potential x'^ arises from a Taylor expansion 
near the bottom of general smooth well. In our work, we have to add a small 
linear perturbation V to the harmonic potential in order to avoid resonances 
(see the non resonance condition ()2.3p below). 

The generalisation of such a result in a multidimensional setting is not 
evident for a spectral reason: the spectrum of the linear part is no more 
well separated. We could expect to adapt the tools introduced in [6] but the 
arithmetic properties of the corresponding spectra are not the same: in |6] 
the free frequencies are Ji + + ' ' ' + Jd jir " ^JD £ while in our 

case they are 2(ji + j2 + • • • + jo) + D for all j'l, • • • , G N. Nevertheless we 
mention that it is still possible to obtain a Birkhoff normal form for (|1.5|) as 
recently proved in [U]. 

A consequence of Theorem 11.11 is the existence of periodic solutions to p.7p . 
There are other approaches to construct periodic solutions of this equation. 
For instance, the gain of compacity yielded by the confining potential allows 
the use of variational methods. We develop this point of view in the appendix. 

The second application concerns the reducibility of a linear harmonic oscil- 
lator, T = -dl + x^, on L^(M) perturbed by a quasi periodic in time potential. 
Such kind of reducibility result for PDE using KAM machinery was first ob- 
tained by Bambusi & Graffi [Tj for Schrodinger equation with an x'^ potential, 
/3 being strictly larger than 2 (notice that in that case the exponent d > 1 in 
the asymptotic of the frequencies (|1.3p ). This result was recently extended by 
Liu and Yuan [14j to include the Dufhng oscillator. 

Here we follow the more recent approach developed by Eliasson &: Kuksin (see 
[7]) for the Schrodinger equation on the multidimensional torus. Namely we 
consider the linear equation 

idtu = —d'^u + x^u + eV{toj^ x)u, u = u{t, x), x G IR, 

where e > is a small parameter and the frequency vector uj of forced oscil- 
lations is regarded as a parameter in U C M". We assume that the potential 
V : T" X ]R 9 (0, x) I— )■ M is analytic in 6 on |Im^| < s for some s > 0, and 
in X, and we suppose that there exists 6 > and C > so that for all 
e G [0, 27r)" and X G M 

(1.8) \v{e,x)\<c{i + x^r^, \d,v{e,x)\<c, \d,.v{e,x)\ <c. 
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In Section [7] we consider the previous equation as a linear non-autonomous 
equation in the complex Hilbert space L2(M) and we prove (see Theorem 17.11 
for a more precise statement) 

Theorem 1.2. — Assume that V satisfies (|1.8|) . Then there exists eo such 
that for all < e < eo there exists Ag C [0, 27r)" of positive measure and 
asymptotically full measure: Meas^Agr) — )• (27r)"' as e — >• 0, such that for all 
oj £ A^, the linear Schrodinger equation 

(1.9) idtu = —d'^u + x'^u + eV{tuj,x)u 

reduces, in L'^(R), to a linear equation with constant coefficients (with respect 
to the time variable). 

In particular, we prove the following result concerning the solutions of (jl.9p . 

Corollary 1.3. — Assume that V is C°° in x with all its derivatives hounded 
and satisfying (jl.Sp . Let p > and uq £ T-L^. Then there exists eq > so that 
for all < e < Sq and a; G A^, there exists a unique solution u£C{R;'HP) of 
(|1.9p so that n(0) = uq. Moreover, u is almost-periodic in time and we have 
the hounds 

(1 - eC)\\u4Hv < \Ht)\\Hv < (1 + eC)\\uQ\\uv. G M, 
for some C = C{p,uj). 

Remark I.4. — In the very particular case where V satisfies (jl.Sp and is 
independent of 9, the result of Corollarv 11.31 is easy to prove. In that case, the 
solution of (II. 9p reads 

u{t,x) = ^c„e"*^"Vn(2;), 

ra>0 

where {(pn)n>o stnd (A„)„>o are the eigenfunctions and the eigenvalues of 
—d'^ + + eV{x), and some (c„)n>o G C. The result follows thanks to the 
asymptotics of ipj when e — )• (see Section [6] for similar considerations.) 

The previous results show that all solutions to (|1.9p remain bounded in 
time, for a large set of parameters u G [0, 27r)". A natural question is whether 
we can find a real valued potential V, quasi-periodic in time and a solution 
u G Ti^ so that ||u(t)||^p does not remain bounded when t — > +00. J.-M. 
Delort [3] has recently shown that this is the case if V is replaced by a pseudo 
differential operator : he proves that there exist smooth solutions so that for 
all p > and t > 0, ||ti(t)||-^p > ct^/^, which is the optimal growth. We also 
refer to the introduction of ^ for a survey on the problem of Sobolev growth 
for the linear Schrodinger equation. 
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Another way to understand the result of Theorem 11.21 is in term of Floquet 
operator (see [5] and |23j for mathematical considerations, and ^ I21j for the 
physical meaning). Consider on L^(M) L^(T") the Floquet Hamiltonian 



then we have 

Corollary 1.5. — Assume that V satisfies (jl.Sp . There exists Eq > so that 
for all < e < Eq and u S A^, the spectrum of the Floquet operator K is pure 
point. 

A similar result, using a different KAM strategy, was obtained by W.M. 
Wang in [23j in the case where 



where hi is the first Hermite function. 

At the end of Section [7] we make explicit computations in the case of a 
potential which is independent of the space variable. This example shows that 
one can not avoid to restrict the choice of parameters to a Cantor type set 
in Theorem 11.21 

Acknowledgements. — The first author thanks Hakan Eliasson and Serguei 
Kuksin for helpful suggestions at the principle of this work. Both authors thank 
Didier Robert for many clarifications in spectral theory. 



We give in this section our abstract KAM result. 
2.1. The assumptions on the Hamiltonian and its perturbation. — 

Let n G M" be a bounded closed set so that Meas(n) > 0, where Meas denote 
the Lebesgue measure in R". The set 11 is the space of the external parameters 
^. Denote by A^^ the difference operator in the variable ^ : 



(1.10) 




n 




k=l 



2. Statement of the abstract result 



Ac,/ = /(-,e) -/(•,??)• 



KAM FOR THE QUANTUM HARMONIC OSCILLATOR 



7 



For I = (Zi, . . . , /fc, . . . ) G so that only a finite number of coordinates are 

oo oo 

non zero, we denote by |/| = \lj\ its length, and (Z) = 1 + | jZj|- We set 

i=i i=i 
Z = {{k, I) / 0, |/| < 2} C X 

The first two assumptions we make, concern the frequencies of the Hamil- 
tonian in normal form (11.211 



Assumption 1 (Nondegeneracy) . — Denote by u = (cji, . . . , a;„) the in- 
ternal frequencies. We assume that the map i— )• a;(.^) is an homeomorphism 
from n to its image which is Lipschitz continuous and its inverse also. 
Moreover we assume that for all (k, I) €z Z 

(2.1) Meas(^{ i : k ■ a;(0 + / • ^H) = }) = 0, 
and for all ^ 

rj](C)/o, VI < 1^1 < 2. 

Assumption 2 (Spectral asymptotics). — Set J^o = 0. We assume that 
there exists m > so that for all i,j>0 and uniformly on H 

— ^j\ > rn\i — j\. 

Moreover we assume that there exists /3 > such that the functions 

are uniformly Lipschitz on U for j > 1. 

If the previous assumptions are satisfied (and actually without assuming 
(|2.ip ). J. Poschel |16j proves that there exist a finite set X C Z and 11^ C 11 
with Meas(n\nQ) — > when a — > 0, such that for all ^ G Ha 

(2.2) \k.u{0 + i-m\>aYTW' {k,i)ez\x, 

for some large r depending on n and /?. 

Then assuming (|2.ip . J. Poschel proves [16] Corollary C and its proof] that 
the non resonance condition (12.21) remains valid on all Z, i.e 



(2.3) \k.Lo{0 + l-n{0\>a-^^, {k,l)eZ,^GUa. 

' ' l + \ky 

In the sequel, we will use the distance 

\^ - n'\2p,U = sup SUp//3 m.(^) _ 0;.(e)| 

fen j>i 
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and the semi-norm 



Finally, we set 



\n\^f,^u = sup supff" \^vM 
C'jen j>i I? - V\ 



kin + l^^lM,n = 



wiiere = sup max — - 



Remark 2.1. — The proof of ()2.3p crucially uses the control of the Lipschitz 
semi-norm |f^l2/3 n (^^^ |161 Lemma 5]). For this reason in assumptions 3 and 4 
below we have to control the Lipschitz version of each semi-norms introduced 
on P or Xp. 

Recall that the phase space V is defined by with a weight ^' so that 

^(j) ^ j; as in the beginning of the introduction. As in [16], for s,r > we 
define the (complex) neig hbourhood of T" x {0,0,0} in V. 

(2.4) D{s,r) = {{9,y,u,v) e V s.t. \lm9\ < s, \y\ < r^, ||n||,i, + \\v\\<i, < r}. 
Let r > 0. Then for W = {X, Y, U, V) we define 

\w\r = \x\ + hri + -( 11^711* + \\vh). 

The next assumption concerns the regularity of the vector field associated to 
P . Denote by 

Xp = {dyP, -deP, d,P, -duP). 

Then 

Assumption 3 (Regularity). — We assume that there exist s,r > so 
that 

Xp : D{s,r) xU^V. 

Moreover we assume that for all ^ G H, Xp{-,^) is analytic in D(s,r) and that 
for all w G D{s,r), P{w, •) and Xp{w, •) are Lipschitz continuous on U. 



We then define the norms 



and 



\P\\D{s,r) ■= sup |P| < +00, 
D{s,r)xn 



, ,,r I ^fn P 

\P\\Dis,r) = sup sup 



n^n D{s,r) I? - 'HI 
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where A^^P = P{-,£,) — P{-,rf) and we define the semi-norms 

\\Xp\\r,D{s,r) ■= sup \Xp\r < +00, 
D(s,r)xn 

and 

\\^P\\r,Dis,r) sup SUp — — < +00. 

C,»)Gn D(s,r) I? - V\ 

where A^^Xp = Xp{-,e) - Xp{-,r]). 

In the sequel, we will often work in the complex coordinates 

1 . ■ N - 1 / 
z = —={u — iv), z = —=(u + iv]. 

Notice that this is not a canonical change of variables and in the variables 
{9, y, z,z)£V the symplectic structure reads 



d9j A dyj + i dzj A 



■J' 



and the Hamiltonian in normal form is 



(2.5) 



i=i j>i 



As we mentioned previously we need some decay on the derivatives of P. We 
first introduce the space '^rD{sry /? > 0, we say that P G ^rD{sr) 

{P)r,D{s,r) + {P)r,D(s,r) < °° "^^^^^ " 

• The norm ( • )r,D{s,r) is defined by the conditions 

2 



(1) 



max 

i<i<" 



If 

dP 



D{s,r) - ^ iP)rMs,r), 



dyj 



D{s,r) 



dP 



dwj 



d^p 



dwjdwi 



< Tg(-P)r,D(s,r), Vi > 1 and Wj=Zj,Zj, 
D{s,r) ]P 

^(^^^ < ■^(P)r,D(s,r), Vj,/>1 and Wj = Zj,Zj. 



'^'Xhis means that {P)r,D(s,r) is the smallest real number which satisfies the mentioned 
conditions : this defines a norm. 
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The semi-norm ( • is defined by the conditions 



max 

i<j<™ 



dP 



dyj 



D{s,r) 



D{s,r) 



^ (-f)r,D(s,r)' 
(s,r)' 



dP 



dwj 

d^p 



dwjdwi 



D{s,r) 
C 

D{s,r) 



< 



< 



Vj > 1 and Wj 



{P)rMs,r)^ Vj,Z>l and 



Zj, Zj. 



The last assumption is then the following 



Assumption 4 (Decay). 



P ^ r r-,/ ^ for some B > 0. 

r,D(s,r) •' ' 



Remark 2. 2. — The control of the second derivative is the most important 
condition. The other ones are imposed so that we are able to recover the last 
one after the KAM iteration (see Lemma l3.4p . Furthermore the assumptions 
on the first derivatives are already contained in Assumption [3] as soon asp > 0. 

2.2. Statement of the abstract KAM Theorem. — 



Recah that M = \uj\^ + \n\ 



2/3,n- 



Theorem 2.3. — Suppose that N is a family of Hamiltonians of the form 
()2.5p on the phase space V depending on parameters ^ G 11 so that Assump- 
tions [i] and are satisfied. Then there exist eo > and s > Q so that every 
perturbation H = N + P of N which satisfies Assumptions [21 and [7] and the 
smallness condition 

e = (||Xp||r,D{s,r) + (-P)r,D(s,r)) + p\\r,D(s,r) + iP)r,D{s,r)) ^ ^O", 

for some r > and < a < 1, the following holds. There exist 

(i) a Cantor set IIq, C 11 with Meas{Il\Ila) —^0 as a — t- ; 

(ii) a Lipschitz family of real analytic, symplectic coordinate transformations 
<^ : D{s/2,r/2) x H^, ^ D{s,r) ; 

(Hi) a Lipschitz family of new normal forms 

n 

3=1 i>i 
defined on D{s/2,r/2) x IIq, ; 
such that 

Ho<^ = N* + R* 
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where R* is analytic on D{s/2,r /2) and globally of order 5 at T" x {0,0,0}. 
That is the Taylor expansion of R* only contains monomials y^z'^z^ with 
2\m\ + \q + q\ > 3. 

Moreover each symplectic coordinate transformation is close to the identity 

(2.6) \\^ - Id\\r,Dis/2,r/2) <Ce, 

the new frequencies are close to the original ones 

(2.7) \oj* -uj\u^ + \n* -n\2p^u^ <C£, 
and the new frequencies satisfy a non resonance condition 

(2.8) \k . oj^O + I ■ > f Y+W' ^) ^ ^' ^ ^ "a. 

As the consequence, for each ^ G Ha the torus $(T" x {0,0,0}) is stUl 
invariant under the flow of the perturbed Hamiltonian H = N + P, the flow is 
hnear ( in the new variables) on these tori and furthermore all these tori are 
linearly stable. 

2.3. General strategy. — 

The general strategy is the classical one used for instance in |10|, I16j . For 
convenience of the reader we recall it. Let H = N+P be a Hamiltonian, where 
is given by (j2.5p and P a perturbation which satisfles the assumptions of 
the previous section. We then consider the second order Taylor approximation 
of P which is 

(2.9) R= Yl Y^Rkm^ije^'-'y^z'^z'^, 

2|m| + |g+g|<2 k€Z" 

with Rkmqq = Pkmqq and wc dcfluc its mean value by 

[R] = Yl Ro^qqy"'^'^- 

\m\ + \q\ = l 

Recall that in this setting z,'z have homogeneity 1, whereas y has homogeneity 
2. 

Let F be a function of the form ()2.9p and denote by Xp the flow at time t 
associated to the vector field of F. We can then define a new Hamiltonian by 
HoXp := A^_|_+P+, and the Hamiltonian structure is preserved, because Xp is 
a symplectic transformation. The idea of the KAM step is to find, iteratively, 
an adequate function F so that the new error term has a small quadratic part. 



12 



BENOIT GREBERT & LAURENT THOMANN 



Namely, thanks to the Taylor formula we can write 

HoX}, = N o Xl + [P - R) o Xl + Ro X], 

= N + {N,F}+ [ il-t){{N,F},F} oXj,dt + 
Jo 



-{P - R)oX}, + R+ [ {R,F}oX^pdt. 
Jo 



In view of the previous equation, we define the new normal form by N-^ = 
N + N, where satisfies the so-called homological equation (the unknown 
are F and A^) 

(2.10) {F,N} + N = R. 

The new normal form A^_|_ has the form (|2.5p with new frequencies given by 

6^+(0 = ^(0 + ^iO and = + 

where 

(2.11) D,(C) = 1^(0,0,0,0,0) and %(e) = ^£-(0,0,0,0,0- 

Once the homological equation is solved, we define the new perturbation term 
P+ by 

(2.12) P+ = {P-R)oX^p+ [\ R{t),F } o X'p dt, 

Jo 

where R{t) = (1 — t)N + tR in such a way that 

HoX}, = N+ + P+ . 

Notice that if P was initially of size e, then R and F are of size e, and the 
quadratic part of P+ is formally of size e"^. That is, the formal iterative 
scheme is exponentially convergent. 

Without any smoothing effect on the regularity, there is no decreasing prop- 
erty in the correction term added to the external frequencies ()2.1ip . In that 
case it would be impossible to control the small divisors (see ()2.3p ) at the next 
step. In this work the smoothing condition ()1.4p on Xp is replaced by As- 
sumption m (see also Remark 14. 3p . The difficulty is to verify the conservation 
of this assumption at each step. 
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Plan of the proof of Theorem 12.31 — In Section [3] we solve the homo- 
logical equation and give estimates on the solutions. Then we study precisely 
the flow map Xp and the composition H o Xp . In Section H] we estimate the 
new error term and the new frequencies after the KAM step, and Section [5] 
is devoted to the convergence of the KAM method and the proof of Theorem 

ESI 

Notations. — In this paper c, C denote constants the value of which may 
change from line to line. These constants will always be universal, or depend 
on the fixed quantities n,/3,n,p. 

We denote by N the set of the non negative integers, and N* = N\{0}. For 

oo 

I = (Zi, . . . , /fc, . . . ) € Z°°, we denote by \l\ = \ lj\ its length (if it is finite), 

oo 

and (1) = I + I'^jljl- We define the space Z = {{k,l) ^ 0, k e Z",/ e 

|/| < 2 }. The notation Meas stands for the Lebesgue measure in M". 

In the sequel, we will state without proof some intermediate results of |16j 
which still hold under our conditions ; hence the reader should refer to |16j for 
the details. For the convenience of the reader we decided to remain as close 
as possible to the notations of J. Poschel. 



3. The linear step 

In this section, we solve equation (|2.10p and study the Lie transform Xp. 
Following [16j, || • ||* (respectively (•)*) stands either for || • || or || • ||^ 
(respectively ( • ) or ( • )^ ) and || • H"^ stands for || • || + A|| • H"^. 

3.1. The homological equation. — 

The following result shows that it is possible to solve equation ()2.10p under 
the Diophantine condition ()2.3p . 

Lemma 3.1 ([16]). — Assume that the frequencies satisfy, uniformly on Ua, 
for some a > the condition (|2.3p . Then the homological equation ()2.10p has 
a solution F, N which is normalised by [F] = 0, [A^] = N , and satisfies for 
alio <a < s, andO< X< a/M 

ll^ivll*,D{s,r) ^ ll^^?llr,D{s,r)' II llr,D(s-<T,r) ^ II^H |lr,D(s,r) ' 

where t only depends on n and r. 
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The space T^^^^^^ is not stable under the Poisson bracket. Therefore 



we need to introduce the space F p^__^ C ^^^^>/.„^ endowed with the norm 



" r,D{s,r) ^ ^ r,D{s,r) 

( ' )rD(s r) + ( ■ ^tois r) defined by the folbwing conditions. 



\D(s,r) — \ /r,D(s,r)' 



max 



dF 



D{s,r) 



<{F) 



r,D{s,r) ' 



dF 



d^F 



dwjdwi 



D{s,r) 



D{s,r) 



< 



< 



j/3+l \ 'r,D{s,r) 
1 



5 Vj > 1 and Wj = Zj, Zj, 



{jini + \j-i\) 



This definition is motivated by the fohowing result, which can be understood 
as a smoothing property of the homological equation 

Lemma 3.2. — Assume that the frequencies satisfy (|'2.3p . uniformly on IIq,. 
Let F,N be given by Lemma \3.1[ Assume moreover that R G '^rD{sr)' ^^^'^ 
there exists C > so that for any < a < s, we have F € ^^'ois-a r)' 
N eV^nf ^ and 

r,lJ(s—cr,r) 



r,D{s—a,r) 



< 



c 



{R)r, 



(3.1) 
and 



C 



r,D{s—(T,r) 



M 



a 



D{s,r) 



(iV)r,D(s-<7,r-) < {R)r,D{s,r), WrM^-^.r) ^ (^)nZ)(s,r)' 

where t only depends on n and r. 

For the proof of this result, we need the classical lemma 

Lemma 3.3. — Let f : M — > C be a periodic function and assume that f 
is holomorphic in the domain \Im9\ < s, and continuous on \ImO\ < s. Then 
there exists C > so that its Fourier coefficients satisfy 

\m\<Ce~\''\' sup \fie)\. 

\Im8\<s 



Proof of Lemma \3.<A — In [16| , the author looks for a solution F of (|2.1U|) of 
the form of (|2.9|) . i.e. 



(3.2) 



F — ^ ^ Fkmqq 

2|m| + |q+g|<2 fcSZ" 
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A direct computation then shows that the coefficients in (j3.2p are given by 

if \k\ + \q-q\^Q, 



(3.3) iFi 



kmqq 



Rkmqq 



k ■ oj + {q — q) ■ Vt^ 
0, otherwise, 



and that we can set = [R\. 

In the following we will use the notation qj = (0, • • • , 0, 1, 0, • • • ), where the 1 
is at the j^^ position, and qji = qj + qi. 

The variables z and z exactly play the same role, therefore it is enough to 
study the derivatives in the variable z. 

In the sequel we write = l + \kY . Then it easy to check that for any j > 1 
and (7 > 0, 

for some C > and t = 2jT + n + 1. In the sequel, t may vary from line to 
line, but will remain independent of a. 

C 

4 We first prove that {F)^^j^(^s_„^r) < — j{R)r,D{^,r)- 

• Observe that 7- — - — = >^ -RfeOa < e**^'^, then according to Lemma [3.31 
dzjOzi ^ ^ 

j.)e~''^'* 

there exists C > so that l-RfcOg^jol < C ' ' , and thus by ()3.3p 

U'J 

and (ESD 



Therefore, as we also have 

(3.5) = ft„,^„e-" 

■> kez^ 

we deduce that 

< 2^ \^kOq,io\e' 



dzjdzi 



D(s—a,r) 



(3.6) < ^(^)rMs,r) 



aat{jiy{l + \j-l\)- 
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• We compute 

(3-7) 1^ = E ^^o,, oe^'-^ + Yl F^^mA'-^i + 2 E ^^'0 2,, oe'^'^z,. 



A>1 



,dR 



Now observe that — )|2^j^q = RkOq^o^^^^ , then by Lemma [33 



|iifeOg,o| < Ce-I'^l^ sup 

|Ime|<s 



< Ce-l^l' 



9i2 



D(s,r) jP ' \ ' / 



From the previous estimate, p.3p and ()2.3p we get 

Ak Cr^fce"!^!'^ 



and thus 



D{s—a,r) 



< 



Jfc|(s-,7) 



(3.8) 



< 



Cr{R)r,D{s,r) 



Similarly, we have |-Ffc02g ol ^ wi r{R)r D(s r)^ which leads to 



(3.9) 



E^^ 



2gj0e 



ik-i 



D{s-a,r) ~ aa*jl^{l + j)' 



< 



By Cauchy-Schwarz in the variable I and 



(3.10) 



D{s—a,r) 



< (E^^'wiE^'^ 



Jk e\2 \ 2 



< 



< 



i>i 

Cr 



fcGZ" 



«>1 



(E /2/3^2(^)(l + _ /|)2) (^)r,D(.,r) 



1>1 

Cr{R)r^D(s,r) 
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since ^(/) > 

Finally, inserting ([331), and (|3TTI]l in we obtain 



(3.11) 



dF ^ Cr{R)r,D(s,r) 



dzj 



D{s-cr,r) aa^j^{l+j) 



dF 



We can write — = V Fkm ooe''''^- Hence by (Q and (12:3 



A 



|-Pfcm,oo| < — l-RfcrnjOol) a^d thanks to Lemma 13.31 applied to the series 
OR 



ik-i 



k€Z" 



(3.12) 



_D(s,r)) 



and we obtain 



(3.13) 



dF 



, < -^fcm,00 e < T{R)r,D(s,r)- 

D(s-a,r) ^-^ ^ aa^ ' V > y 



To obtain the bound for ||-F'||D(s-(T,r) write 



(3.14) F= 5^ Ffcoooe''^-^+ ^fc™. ooe'*^-% + 

fceZ",i<i<n 



jk-e- 



k&Z",j,l>l 



k€Z",j,l>l 



k&Z",j,l>l 



Since R^y^^^-^^Q = Y^ -Rfcoooe'^^) by Lemmas 13.31 and 13.11 we deduce that 



(3.15) 



|F.ooo| < Cr'^e-\''\^R) 



r,D{s,r) ) 



hence, thanks to ()3.12p and ()3.15p we can bound the sums of the first line in 
()3.14p as in the previous point. 
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Now thanks to ()3.4p and to the Cauchy-Schwarz inequaUty we have 

/ ^ ^kOq.iOe ZjZi < > 



< 



D(s-a,r) OCT* + \j - l\ 



C{R)r,D{s,r) f s---^ \Zj\^'^ 



Cr^{R)r D(s r) 

Therefore we proved that iDfs-o-r") < 7 ' — • 

This latter estimate together with the estimates (j3.6p , (j3.1ip and (jS.lSp shows 
that 

C 

^ We now show that 

(3.16) {N)r,D{s~a,r) < {R)r,D{s,r)- 

Since N = [R] we have 

n 

(3.17) N = ^Ro^ 

and we can observe that 
(3.18) 

1 f OR 1 f O^R 

= ^(«,o.o,o)d«. = ^ _(«,o,o,o)de. 

which im ply th e bounds |-RomjOo| < {R)r,D{s,r) and l-Roog^g,! < {R)r,D{s,r)/f^ 
and thus ()3.16p . 

^ It remains to check the estimates with the Lipschitz semi norms. 

As in |16j . for \k\ + \qj — ^| 7^ define 6kji = k ■ oj + il.j — Qi. Then by ()3.3p . 

iA^r]Fkmqjqi = ^kjl^^^^^rjRkmqjqi + RkmqjqiiO^^ri^kJl- 

By (I23D, 1(5;:]; I < Ak/a and thus 

^2 
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hence 



1^ — r]\ ' ' 

and we have 

(3-19) |e-r?| -^^l |e-r?| ^ ^\^^^'^f^M)- 

Thanks to the estimate ()3.19p it is easy to obtain ()3.ip . 

Finahy, the estimate {^)j.^D{s-cr,r) (s,r) ^ straightforward conse- 

quence of (IXrm and ^M)- □ 

3.2. Estimates on the Poisson bracket. — 



r,U(s,r) r,U(s,r) 

of the form (|2.9p . Then there exists C > so that for any < a < s 



Lemma 3.4- 



Let R e ri;'n,_x and F G r„'^/„_N be both of degree 2, i.e. 



(3.20) 



C 



{{R,F} )r,D{s-a,r) < " (^>r,D(s,r) (-^)^z)(,,r) 



anc 



Proof — The expansion of { i?, F } reads 

.^^.A/aRap aRaF\ -^(^^ dn dF \ 



^ ydOk dyk 

It remains to estimate each term of this expansion and its derivatives. We will 
control the derivative with respect to 9k thanks to the Cauchy formula : 



(3.21) 



dP 



<-\\P\\ 



which explains the loss of a. 

Notice that if P is of degree 2 (and that is the case for F and R) we have 

Q2p Q2p Q3p 



(3.22) 



dzdy dy^ dz^ 



0, 



fact which will be crucially used in the sequel. Finally observe that z and z 

d 

exactly play the same role, hence we will only take — into consideration. 

oz 

4 We first prove (Km . 
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a 



{R)r,D{s,r){F)j^D{s,r)- 



• Since ll-P Q||D(s,r) ^ ll-P||D(s,r) IIQIlD(s,r) we have by Cauchy formula 

k>l 

, , Cr^ 
(3.23) < 

• With (fOTT) we have 



D{s,r) 



c 

estimate holds interchanging R and F. In view of (|3.22|) 





fdRdF\ 


< 

D{s-a,r) 




fdR\ 




dF 




KdOk dyk) 


dOk 


Kdyj) 


D{s—cr,r) 


dyk 



C 


dR 




dF 




a 


dyj 


D{s,r) 


dyk 


D{s,r) 



and the same 
deduce 



we 



d C 
(3.24) max —{R,F] < -{R)r D(s r){F)^ n, ^ 

^ ' l<y<n dyj^ ' J Dis,r)- /r,V(s,r)\ /r,D{s,r) 

• By ([3:22]) ^ /dR dF\ _ OR d^F 

' dzjKdykdOk) dykdzjdO^' 



dR 


d^F 




dyk 


dzjdOk 


D{s-a,r) 



c 


dR 




dF 




a 


dyk 


D{s,r) 


dzj 


Dis,r) 



< 



d (dRdF 
Similarly ^(^^ 
dzj V d0k dyk 

rule 

d rdRdF\ 

dZj \ dZk &Zk / D(s,r) 



< 



Cr 



'\R)r,D{s,r){F)XD(s,ry 



Cr 

D[s-a,r) - -jF^^^^-'D(^'-)^^^rMs,rr t^C Leibniz 



< 



d^R 




dF 


+ 

D{s,r) 


d^F 




dR 


dz^dzj 


D{s,r) 


dzk 


dzjdTzk 


D(s,r) 


dzk 



Cr /I 1 \ 

- i^VA^ ^ fc2/3(l + \j - fc|) j^^^'^'^(«''^)^^^^0(s,r)' 

and taking the sum in k yields 

^\\±(dRd^\ <CV 

^ WdZj \dZk dzk) Dis,r) - j^a^^'"'''^'''-^^ '^'D^^'^y 



D{s,r) 



k>l 
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The previous estimates imply that 
d 



(3.25) 



dzj 



{R,F} 



r,D{s,r){F)r,D{s,r)- 



^, , , , 92 fdRdF\ dR d^F „ , 

Thanks to (^22), 7r^\l^7^) = ^ ' (^ZB) we 

dzjdzi \ dyk oOk^ dyk dzjdzidOk 

(- 



obtain 



/ dR OF 



dz-jdzi \dyk 80 u 



D(s—a,r) 



(3.26) 



< 



< 



dR 




d^F 




dyk 


D{s,r) 


dzjdziddk 


D{s~a,r) 



c 



{R)r.D(s,r){F) 



,D(s,r)\^ lr,D(s,r)' 



and the same estimate holds interchanging R and F 
On the other hand, 

d"^ / dR dF 



d'^R d'^F d'^R d'^F 
+ 



dzjdzi \dzk dzk' dzjdzk dzidzk dzidzk dzjdzk^ 



and 



d^R d^F 



dzjdzk dzidzk 



D(s—(7,r) 



< 



< 



d^R 




d'^F 




dzjdzk 


D{s,r) 


dzi&Zk 


D{s,r) 



c 



{jlk^Yil + \l- k\) ^^^^'^(«''^)^^^r,D(s,r)- 



Hence, with (|3.26p we conclude that 
(3.27) 

II UZjUZl ' 

1 



q2 (J 



as the series 



fc>i 



k'^l^{l + \l-k\) 



converges. 



Finally, the estimates ([323]), (lOiD . ([3:25]) and IKTi}i yield the estimate IKlQh . 

^ To prove the estimate with the Lipschitz norms, we can use the previous 

analysis and the two following facts. 

Firstly, since A^rjifg) = f{0^Cr,9 + g{r])Ai:,^f, hence 

\\f9\\D(s,r) ^ \\f\\D{s,r)\\9\\Dis,r) + ll^ll D{s,r) II / ll£(s,r) • 

Secondly, the operator A^^ commutes with the derivative in any variable. □ 
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3.3. The canonical transform. 



In this Section we study the Hamiltonian flow generated by a function F S 



" r,D{s—<T,r) 



globally of degree 2, i.e. of degree 2 in the variables z, z and of 



degree 1 in the variable y. Namely, we consider the system 

r (^(i),y(t),i(t),i(0) =X^((e(t),y(t),z(t),z(t))), 
I (e(0),y(0),z(0),z(0)) = (00,^0,^0,^0). 



(3.28) 



Lemma 3.5. — Let < a < s/3 and F G '^r'D{s-cTr) ^^^^ ^ degree 2. 
Assume that {F)^jy^^_^^^ < Ca. Then the solution of the equation ()3.28p with 

initial condition (^9^ , , jZ^) G D{s — 3a,j), satisfies (^6{t),y{t), z{t),^{t)^ G 
D{s — 2a, |) for all < t < 1, and we have the estimates 



(3.29) 
(3.30) 



sup 

0<t<l 



sup 

0<t<l 



dyk{t) 



< 



Cr{F) 



r,D{s—a,r) 



with Wj 



z] or^^, 



dwk{t) 



dwj 



< 



\ lr,D{s-a,r) , r -.r _ „ 

(jfc)/3(l + |j-A:|) ^ ^'^ " '''' ""^ " '"^ 



(3.31) 
(3.32) 



sup 

0<t<l 



sup 

0<i<l 



dyk{t) 



dy^ 



< 



D{s-a,r) 



a 



dhkit) 



< 



C{F) 



r,D{s—a,r) 



a{ij)('{l + \i-j\) 



with w\ 







zO orz?, upj 



zO or zO. 



Before we turn to the proof of Lemma 13.51 we introduce a space of infinite 
dimensional matrices, with decaying coefficients. 

Let II • II be any submultiplicative norm on Al2,2(C), the space of the 2x2 
complex matrices. For /3 > 0, we say that B G A^f if (( B ))^^ < oo, where 
the norm (( • is given by the conditior^ 



sup sup ||i?j7|| < 



{jif{^ + \3-i\y 

Then we have the following result 



Vj,/ > L 



'^^This means that {( • ))/3 s smallest real number which satisfies the mentioned condi- 

tions : this defines a norm. 
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Lemma 3.6. — Let A,B e M 



Then AB £ M 



0, 



and 



k>l 



Proof. — For all j,l > 1, [AB)^^ = ""^^ Aj^Bki. Since || • || is submultiplicative 
\\{AB).^\\ < Y^WAjkWWBkiW 

k>l 



(3.33) 



< 



Thanks to the triangle inequality, for all j,l > 1, 

{A:>1}C{A:>1 : | j - fe| > - /| } U > 1 : \l - k\ > - l\} , 
thus, by splitting the sum in (|3.33|) we obtain the desired result. 



Proof of Lemma \K 
Z = iZj)j>i. Then F reads 



□ 



Here we introduce the notations Zj = {zj,Zj) and 



1 



(3.34) F{e,y,Z) = bo{0) + h{9)-y + a{e)-Z + -{A{d)Z)-Z, 

with 

boie) = F{e, 0, 0), h{9) = VyF{9, 0, 0), a{e) = VzF{e, 0, 0), 
and A = (Aij) is the infinite matrix so that 

/ d^F 

(3.35) 



dzidzj 



(9,0,0) 



^{o,o,oy 



dzidzj 



Q2p Q2p 

\ dzi dzj ^^^^'^^ Q-^. Qtz- ' ^ ' 



Observe that A is symmetric. 

By [lUl Estimate (9)], the flow exists for < t < f and maps D{s — 3cj, |) 
into D{s — 2a, |). Here we have to give a precise description of X^p for < 
t <1. This is possible thanks to the particular structure (j3.34p of F. 
In the sequel we write {e{t),y{t), Z{t)) = X*,((9°, Z^). 

l|k To begin with, the equation for 9 reads 



(3.36) 



9{t) = VyF{9,0,0)=bi{9), 9{0) 
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Since bi is a smooth function (see (j3.2p ). the n-dimensional system (j3.36p 
admits a unique (smooth) local solution 0{t). By the work of J. Poschel, this 
solution exists until time t = 1, and we have the bound 

(3.37) sup |Im6'(t)| < s - 2cj, 

0<t<l 

(this can here be recovered by the usual bootstrap argument, using the small- 
ness assumption on F). 

^ We now turn to the equation in Z. We have to solve 

(3.38) Z{t) = JVzF{e, y, Z){t), Z(0) = Z^, 
where 

J = diag{(_°^ J)},>1- 
Notice that by |16|, Estimate (9)] we already know that 



(3.39) sup \\Z{t)\\p < 

o<t<i * 



r 



but we need to precise the behavior of Z(t). 

Since 6 = 9{t) is known by the previous step, in view of ()3.34p . equation ()3.38p 
reads 

(3.40) Z{t) = b{t) + B{t) • Z{t), Z(0) = 

where b{t) = Ja{e{t)) and B{t) = JA{9it)). 

We now iterate the integral formulation of the problem 

Z{t) = + /* {b{ti) + Bih) ■ Z(ti))dti, 
Jo 

and formally obtain 

(3.41) Z{t) = 6°°(t) + (1 + B°°{t)) ZO, 
where 



(3.42) ^~W = E/„/„ /„ llB{tj)b{tk)dtk---dt2dtu 

and 



(3.43) B^{t) = y2 / ■■■ I l\B{tj)dtk---dt2dh. 

By (f335]) and (13:3711 . there exists C > so that 

sup \\B{t)\\i2 (2 < C, 



0<i<l * * 
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and thus, for all < t < 1 the series ()3.42p converges and 

,2 < sup \m\\,2 T c'-' f r ■■■ dtfc ■■■dt2 dti 

o<t<i * fr: Jo Jo Jo 



k>l 

wmwa. } . 

0<t<l 



fc>i 
e^ - 

< sup \\b{t)U—- 
0<t<l * 



(3.44) < C sup 



o<t<i * 



Similarly we have uniformly in < t < 1 



||i?°°(t)||,l^,l<a 

As a conclusion, the formula (j3.4ip makes sense. 

Indeed, we need more precise estimates on Recall that B[t) = A[9{t)), 

where A is defined by Then by and (I23Z1), for all < t < 1, 

B{t) G A^fL+ and sup {{B{t)))+ _^ < C{F)+ Hence by Lemma 

0<t<l . V : ; 

[Mland dOHjl 

(3.45) < e^<^>^-(^--) - 1 < 
l|k Finally we turn to the equation in y 

yit) = -VeF{9,y,Z){t), y(0) = y". 

We already know the functions 9{t) and Z{t). Moreover as the function F 
(j3.34p is linear in y, the previous n— dimensional system reads 

(3.46) y{t) = f{t) + g{t)y{t), y(0) = y", 
with 

fit) = -VeboiOit)) + Veaieit)) ■ Z{t) + \[y eA{d{t))Z{t)) ■ Z{t), 

and 

g{t) = -Vebi{e{t)) = -VeVyF{e,0,0). 

We can solve the equation ()3.46p with the same techniques as the equation 
()3.38p . In fact we have formally 

(3.47) y(^t) = rit) + {l + g'^it))y', 
where 

(3.48) rit) = y] f r ■■■ f'~'T\g{t,)f{tk)dtk---dt2dt^, 
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and 



t rt 



Yl9itj)dtk---dt2 dh. 



By (|3.37p and the Cauchy formula 
sup < — max 



dF 



D{s—a,r) 



< 



C{F] 



r,D{s—cr,r) 



a 



and similarly to (|3.44p we have for all < t < 1 



\r{t)\<c sup 

0<t<l 



and 



(3.49) 



< 



D{s—cr,r) 



a 



which shows the convergence of the series defining ()3.47p . 

4k It remains to show the estimates on the solutions of ()3.28p . 

• First we prove ([330]) . By (I33Q|), 

therefore by (j3.45p , for k ^ j we have 

\ ' j',U(s—a,r) 



(3.50) 



and IIV70Z 



< h 



■■■■ ijk)^{l + \j-k\)' 
which was the claim. 

• We prove (jOT]) . By ([HliTP we have 

l<j<n 

hence — -q = 5jk + djki^) the claim follows from ()3.49p (/°° does not 
depend on y^). 

• We prove (|3.29p . Since g and do not depend on Z, from (j3.47p we deduce 
dy df^ 



that 



Now by definition ^Mh of f°°, we get that for all < t < 1 

(3.51) 'y^'^ - 'f^^^'^ 



< Vzof^it) <C sup |V^o/(t)|. 

^ 0<t<l ^ 
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For all 1 < / < n, we compute 

(3.52) ^zjiit) = deM^it)) + ^ de,Aki{e{t))Zi{t). 



i>l 



As ak{0) = V Zf,F{9,Q,Q), with the Cauchy formula we deduce 



sup \de,akm))\ < < 



Cr{F) 



r,D{s—a,r) 



0<t<l 

Similarly with (j3.35p . 



sup \de,Aki{e{t))\ < 



C{F) 



r,D{s—a,r) 



0<t<l 



a{ik)l^{l + \i-k\)' 



Inserting the two previous estimates in ()3.52p . we obtain using ()3.39p and the 
Cauchy-Schwarz inequality 



\^zji{t)\ < 



C {F)t 



D{s—cr,r) 



(3.53) 



< 



Cr{F) 



kP 



r,D{s—cr,r) 



'^iP{l + \k-i\ 



a 



kP 



Since Vzofi{t) = ^ {V zoZk{t))V zji{t), from (1330]) and ([333]) we deduce 

k>l 

iv^o/Koi < Y.\\^zoZkm\\vzjim 

;D{s—a,r) ( ■ST'^ 1 



< 



k>l 

Cr{F)^ 



k> 



.A;2/3(l + |j-A;| 



+ 1 



< 



Cr{F) 



r,D{s—a,r) 



0-J 



/3 



and together with ()3.5ip . we get that for all j > 1 



sup 

0<i<l 



dy{t) 



< 



Cr{F) 



r,D{s—a,r) 



It remains to show ()3.32p . first we have 
dy{t) 



dzfdz^ 

^ J 



< 



< C sup 

0<t<l 



VzoV2o/(t) 
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Then from the very definition of /, V zoVzof{t) = VQAij{6(t)), and using the 

* 3 

Cauchy estimate in 6 we get, 



dy{t) 



< 



C{F) 



r,D{s—a,r) 



a{ijr{l + \i-j\y 



which was the claim. 

In the next result, we denote by | • |^ the Lipschitz norm 



□ 



sup 



1/(0-/(^)1 



We have an analogous result to Lemma 13.51 with Lipschitz norms. 

Lemma 3.7. — Under the assumptions of Lemma \3.5\ and the condition 
r,D{s-(T,r) — solution of (|3.28p satisfies moreover 



sup 

0<t<l 



sup 

0<t<l 



dykit) 



dwk{t) 



sup 

0<t<l 



dyk{t) 



sup 

0<t<l 



dhk{t) 



< 



< 



< 



< 



Cr{F) 



r,D{s—cr,r) 



0-J 



with vJj 



C{F) 



r,D{s—iT,r) 



Ukni + \j-k\) 



a 



C{F) 



r,D{s—<T,r) 



a{ij)l^{l + \i-j\) 



z] or^j, 



with Wk = Zk or z^, 



Zj or z^, 



with w. 







Zi or w° 



z] or z]. 



Proof. — We won't detail the proof, since it is tedious and similar to the 
proof of Lemma 13. 5i First we define the space Ais'~^'^ with norm {{■))~^'^ 
similarly to A^s'"*", but with a Lipschitz norm in C Then we have {{AB))+^^ < 
C(((^))+'^((-B))++((S))+'^((yl))+). Then one can follow the proof of Lemma 
13.51 and use that the different norms (say || • ||) which appear satisfy 1 1/(71 1''' < 



C 



9\n- 



□ 



To conclude this section, we state a result which shows that the Lie trans- 
form associated to a quadratic function, is also quadratic. This will be crucial 
in the proof of Theorem 12.31 (see Section [5.2|) . 
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Corollary 3.8. — The symplectic application Xp reads 
9 \ ( K{6) 

y \ ^{ L{e,z) + M{e)z + s{e)y 
z j \ T{e) + u{e)z 

where L{6,Z) is quadratic in Z, M{6) and U{9) are bounded linear operators 
from i\ X i\ into itself and S{9) is a bounded linear map from M" to W^. 

Proof. — The claim follows from the proof of Lemma 13.51 The structm'e of 
Z{1) follows from ()3.40p . while the structm'e of y(l) comes from (j3.47p and 

3.4. Composition estimates. — 

In this section we study the new Hamiltonian obtained after composition with 
the canonical transformation Xp. 

Proposition 3.9. — Let < rj < 1/8 and < a < s, R ^ '^lyr D(s-2a i-qr) '^^^ 
F e r'^'t . with F of degree 2. Assume that (F)+^, . + (F)^'^. . < Ca. 

r,D(s—a,r) ■' ^ > ' r,D{s,r) > ' r,D(s,r) 

Then R o Xl, G c n an-d we have the estimates 

(3.54) {R° Xp)^r,D{s-5a;r]r) ^ C { R) r]r,D(s~2a,'lrir) ^ 

{RoXp )^r,D(s-5a,rir) ^ C {{ R ) r,r,D{s-2a,irir) + { R )^r,D(s-2<T,47?r) ) • 

Proof. — The proof of the first estimate relies on Lemma 13.51 We omit the 
proof of the second, which is similar using the estimates of Lemma [3.71 instead . 
In the sequel, we use the notation {6, y, z, z) = Xp{9^, y^, z^,z^). 
^ Since Xp maps D{s — 3a, j) into D{s — 2a, §), it is clear that 

(3-55) ^F\\D{s-5a,r]r) ^ C (R) rjr,D{s-2a,4rir) ■ 

^ By the Leibniz rule, for all 1 < j < n 



dy'j ^ dyk 



and by ()3.3ip we deduce 



d{RoXj,) 

l-3-t>t)J — p ^^^{^)r,r,D{s-2a,4vr)- 



dy] 



D{s—5cr,'rir) 
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^ For j > 1, the derivative in reads 



d{R o Xj,) 



\dR{X^p)dyk , ^fdR{X},)dzk , dR{X^p)&Zk 



E 



+ 



fc=l ^'^ 1 k>l 
Therefore, thanks to ([3:29]) and ([3:32]) we get 

\d{RoX^p) 



dzk dz^- ^ cfzk dz^- 



(3.57) 



< 



dz^j D{s-5a,rir) 

7^\\dR{X}, 



< 



E 

k=l 



dyk 



D 



\ dyk \ 

{s-5a,r]r) I dZj I 



fc>i 



D(s-5(T,r) 



dzO 



C ST ^ 



< -Tp{R)-qr,D(s-2(7,irir)' 



l|k We now estimate 



D{s—5cr,r]r) 



for i,j > 1. By the Leibniz rule, 



the result will follow from the next estimations. 
• Using the Cauchy estimate in yi and (|3.29|) 



d^RiX^^) dyk dyi 



dykdyi dz^ dz[ 



< 



By iKm 



^ dRjX^p) d^yk 



l<k<n ' J 



By (IM]) 



E 



d''R{X^p) dzk dzi 



k,l> 



dzk{t)dzi dzi dz^j 



C{R)nr,D(s-2<T,4rir) 
Dis-5a,r,r) ~ {ijY 

C {R)r)r,D(s-2a,Arir) 

w ■ 

C{R)jjr,D{s-2(j,4rir) 

w ■ 



D{s—5(7,rjr) 



< 



D{s—5ij,rir) 



< 



Using the Cauchy estimate in z^, ([3.29P and p.30p we get 



E 



d^RiX'p) dyi dzk 



k> 

KKn 



^ dzkdyi dz^ dz*^ 



D(s—5cr,r]r) 



< 



C{R)nr,D{s~2a,4:rir) 
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Finally, (|33i]) follows from (|335]) . (1336]) . (ITO]) and (f338]) . 



□ 



3.5. Approximation estimates. — 

Recall that the notation || • ||* (respectively ( • )* ) stands either for || • || or || • H"^ 
(respectively ( • ) or ( • )^ ). 

First we recall some approximation results \16\ Estimate (7)], which show that 
the second order approximation of P can be controlled by P, and that P — R 
is small when we contract the domain (this contraction is governed by the new 
parameter r]): 

Lemma 3.10 ( jl6| ). — Let P satisfy Assumption\^ and consider its Taylor 
approximation R of the form (|2.9p . Then there exists C > so that for all 
ri > 



\\^R\\r,Dis,r) < C'l|-'^p||*,D(s,r)' ^'^^ P - ^R\\*rjr,D{s,4vr) < C'^ II -'^P II r,Z)(s,r) ' 



Proof. — • We first prove the second estimate. Define the one variable func- 
tion f{t) = P{6,t'^y,tz,tz). Then by the Taylor formula, there exists < to < 
1 so that 




and 



/(l) = /(0) + /'(0) + -/"(0) + -/(3)(to) 



which reads 



P{0,y,z,z) - R{e,y,z,z) 



1 



f^'\to) 



6 




Using the Cauchy estimates in z or in y, we obtain 



P - ^llD(s,4r,r) < Crj {vrf {P)r,D(s,r)- 
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The estimates of the derivatives are obtained by the same method, with the 
adequate choice of the function /. A derivative in z costs r] and a derivative 
in y costs r/^. 

It is then also clear that we have (P - R)^r,D{s,4.rir) - ^'^(^)r,D{s,r)- 

• The inequality (P)*£)(5r) — ^^^Ko^sr) ^® ^ consequence of the previous 
point with tj = 1. □ 



4. The KAM step 

Let be a Hamiltonian in normal form as in (|1.2p , which reads in the variables 

i<j<" i>i 

and suppose that the Assumptions [T] and [2] are satisfied. 

Consider a perturbation P which satisfies Assumptions O and U] for some r,s > 

0. Then chose < r/ < 1/8, < o" < s, and assume that 

(4.1) {P)rMs,r) + \\Xp\\rMs,r) + ^ + \\Xp\\^Ms,r)) < 

where t is given by Lemmas 13. II and 13.21 cq is a large constant depending only 
on n and r (see [161 Estimate (6)].) 

Thus, by Lemmas 13.11 and 13.101 the solution F of the homological equation 
(fZTOD satisfies 

\\XFrr,Dis-.,r) < ^ II^P II^D(s-.,.) < ■ 

Similarly, by Lemmas 13.21 and 13.111 

SO that the hypothesis Lemma 13.51 are fulfilled. 
We use the notations of Section 12. 3[ 

4.1. Estimates on the new error term. — 

We estimate the new error term P+ given by ()2.12p . 
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Lemma — Assume ()4.ip . Then there exists C > (independent of t] 

(X 

and a) so that for all < X < — 

\ + / rir,D{s—5a,r]r) ' II P+ \\'qr,D{s—5(j,'qr) — 

Proof. — By Estimate (13)], we already have 
(4.2) \\XpJ 

r]r,D{s — 5cr,rir) — rurr^n'^ 

i\\Xp\\ 

r,D{s,r) r.D(s,r)- 

It remains to prove a similar estimate for the (, ) norm. 
By Lemmas 13.51 and 13.111 

((^ - R)° XF)^r,D{s-5a,'nT) - ~ ^)^r,D{s-2a,'irir) - ^V{P)r,D{s,r)- 

Then by Lemma 13.51 again 

(/"'{i2(t),F}oX^dt)J,,^(,_5^^^,) < C A{i?(t),F}oX^)J,^^(,_5^^^,)dt 

•/ 

< C{{R{t),F})^^jy(^^_2^,^^j.y 

Since R ^ and F G T^'^, n are both of degree 2 we can apply 

r,D(s,r) r,D(s—c7,rj 

Lemma 13.41 and write 



( I { ^(Oi F]°X^F d*>^r,D(s-5a,,?r) ^ — (-R)^r,D(s,,?r) (^)^r,D(s-<7,r,r) • 
J 

Finally by Lemmas 13.21 and 13.111 

where we used that {■)r]r,D(s,r]r) ^ 'n~'^{')r,D(s,r)- Putting the previous esti- 
mates together, we complete the proof. □ 



4.2. Estimates on the frequencies. — 

We turn to the new frequencies given by ()2.1ip . 
Lemma 4 •2- — There exists K > 10 and q;+ > so that 

\k-oj+{0 + l-^'^{0\>a+^, \k\<K, \l\<2. 

In fact K can be made explicit, it depends on n, r, cq and on all the constants 

C. 
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dN 

Proof. — On the one hand, since u}j(^) = — — (0, 0, 0, 0, 0, by Lemma r3.1UI we 

9yj 

deduce that 



ON 

\uj\n < sup | — I < ll-'^yyllr.Drs.r) < C\\XR\\j.^D(s,r) < C\\Xp\\j.^D(s,r)- 
D{s,r)xU Oy 

On the other hand, = ^ ^ (0,0,0,0,g), thus 

OZjOZj 

(4.3) 

\^\2p,U< sup I < (A^)r,D(s-a,r) < C{R)r,D{s,r) < C'(-P)r,D(s,r), 

D(s,r)xn OZjOZj 

hence by the two previous estimates 

(4.4) |a5|n + |S^|2/3,n < C'(l|-'^p||r,D{s,r) + {P)r,D{s,r)) ■ 

Similarly, for the Lipschitz norms we obtain 

Pin + M2l3,n < C'(ll-'^P|l^;D{s,r) + {P)r,D{s,T))- 

We follow the analysis done in |16j to bound the small divisors and thanks to 

g31) 

\k-'Cj + l-h\ < |A;|(/)(|a}|n + |f^|2/3,n) 

< C\k\{l){\\Xp\l^nis,r) + {P)rMs,r))- 

We now choose a > CqK max ^fc(||^p||r,D(s,r) + {P)r,D(s,r)) where Co is a 

\k\<K ' ' ' ' 

large universal constant, and thanks to the estimate given by the frequencies 
before the iteration we get 

|fc-a;+(e) + /-J^+(OI >«+^, \k\<K, 

with a+ = a — a. It remains to show that a"*" > 0. This is done in [16\ Section 
4], and the proof still holds with the new norms. □ 



Remark 4- 3- — The key point in the previous proof is the estimate ()4.3p . 
which shows that the perturbations of the external frequencies can be con- 
trolled by {P)r,D{s,r)- In the case of a smoothing perturbation P (case p > p 
in (|1.4p ). the norm {•)r^D{s,r) is not needed (more precisely, the decay of the 

derivatives of P is not needed), because we then have |f^|2/3,n < ll^p||r,D(s,r) 
with /3 = (p-p)/2. 
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5. Iteration and convergence 

In this section we are exactly in the setting of [16], and we can make 
the same choice of the parameters in the iteration. We reproduce here the 
argument of J. Poschel. 

5.1. The iterative lemma. — 

Denote Po = P and Nq = N. Then at the i^— th step of the Newton scheme, 
we have a Hamiltonian Hu = Ny + Pu, so that the new error term Pu+i is 
given by the formula (j2.12p and the new normal form N^^i is associated with 
the new frequencies given by (|2.11|) . 

Let ci be twice the maximum of all constants obtained during the KAM step. 
Set ro = r, So = s, ao = a and Mq = M. For u >0 and k = 4/3 set 

a, = ^(1 + 2--), M, = Mo{2-2-n, = 



Cl£u <^v 3 



and 



The initial conditions are chosen in the following way : ctq = so/40 < 1/4 so 
that So > si > • • • > So/2, 

£0 = 7oao<7o and 70 = (cq + 2*+^ci)~^, 

where cq is the constant which appears in (|4.ip . We also define Ky = Kq2'^ 

with K^+^ = l/(ci7o). 

With the notation Di, = D{sy, r^) we have 

Lemma 5.1 (Iterative lemma, |16)). — Suppose that Hy = Ny + Py is 
given on Dy x Uy, where Ny = u]y{^)-y + ^y{C) ■ zz is a normal form satisfying 

\k-UJy{i)+l-Vty{0\>ay^, {k , I) ^ Z , 

on Uy and 

{P)tD. + \\Xp\\tD.<e.. 

Then there exists a Lipschitz family of real analytic symplectic coordinate 
transformations <&jy+i : Di^+i x 11,^ — > Dy and a closed subset 

n,+i = n,\ J ni+\ay+i), 

\k\>K^ 
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ofHu, where 

such that for H^^i = Hy o = A^i^+i + Pu+i, the same assumptions are 

satisfied with v + 1 in place of v. 

We don't give the details of the proof of this result, since it is entirely done 
in |16| : it is of course an induction on G N which essentially relies on the 
results of the Section HI 



5.2. Proof of Theorem [231 — 



The result of Theorem 12 .31 is the convergence of the sequence H^, to a Hamil- 
tonian in normal form, for parameters ^ in a set IIq,, which is the limit of the 
sets liy. 

We again follow the proof of Poschel and we recall the following Lemma 

Lemma 5.2 (Estimates, [16j). — For u >Q, 

1 Ce 

Set Ho = n\Ufc,i'7^fcf and n„ = n,y>iUy. The proof that Meas(n\nc,) — > 
when a — > is done in |16i, Section 5] and we do not repeat it here. 
For > 1 we define the map 

^.'^ = $1 o . . . o : X n^_i Dy_i, 



and thus we have H,^ = H o With the Lemma 15.21 and since^ Cii^yiDi, x 
Hiy = D{s/2) X Hq, we are then able to show, as in |16) . that is a Cauchy 
sequence for the supremum norm on D{s/2) x Hq. Thus it converges uniformly 
on D{s/2)xlla and its limit <I> is real analytic on D{s/2). Further, the estimate 
([22D holds on D{s/2) x n„. 

It remains to prove that $ is indeed defined on D{s/2, r/2) x with the same 
estimate. By Corollary 13.81 all the transforms are linear in y and quadratic 
in z,z and thus the same is true for the transform $ (this fact was also used 
in|17j or [B]). This specific form is stable by composition and thus all the ^'^ 
have this form and in particular they are linear in y and quadratic in z, z. 
Therefore it suffices to verify that the first derivatives with respect to y, z, z 
and the second derivatives with respect to z, z of are uniformly convergent 



'^'here we use the notation D{s/2) = _D(s/2,0). 
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on D{s/2)xlia to conclude that convergences to (actually an extension of 
the previously defined $) uniformly on D(s/2, p) x ITq for any p. In particular, 
for r small enough, 

$ : D{s/2,r/2) xUa-^ D{s,r) 
and <I> still satisfies estimate (|2.6|) . 

So it remains to analyse the convergence of the derivatives. Using Lemma 15.21 
we obtain successively ||-D<I>;^||r„^j.^^£)^ < 2 and then uniformly on D{s/2) x 

and we deduce that uniformly on D{s/2) x IIq, 

So again D^'^ converges uniformly on D{s/2) x IIq,. Similarly we obtain the 
convergence of the second derivatives using the formula 

On the other hand, again using Lemma 15.21 the frequencies functions ooy and 
converge uniformly on IIq, to Lipschitz functions oo* and satisfying ()2.7|) 
and thus ()2.8p in view of Lemma 15.11 
We then deduce that, uniformly on Z)(s/2,r/2) x IIq, 

:= H o - — > Ho^-N*=:R* 

and since for all v the Taylor expansion of Ry contains only monomials y'^z^z^ 
with 2\m\ + I^ + qI > 3 the same property holds true for R* . □ 



6. Application to the nonlinear Schrodinger equation 

Let n > 1 be an integer and i/, e > be two small parameters so that u > 
CqE, where Co > is a constant which will be defined later. Set 11 = [—1, 1]". 
We consider a perturbation of the one dimensional Schrodinger equation with 
harmonic potential 

(6.1) idtu + dlu-x'^u-iyV{^,x)u = e\uf'^u, {t,x)eRxR, 

where m > 1 is an integer and {V{^, ■))^gn ^® family of a real analytic bounded 
potentials with ^(0, ■) = which will be made explicit below. 
Recall that T = —d"^ + denotes the harmonic oscillator. Its eigenfunctions 
are the Hermite functions (/ij)j>i, associated to the eigenvalues (2j — l)j>i. 
Now consider the linear operator A = A{v,^) = —d'^ + + i'V{^,x). Under 
the previous assumptions, A is self-adjoint and has pure point spectrum with 
simple eigenvalues {Xj{i^,^))j>i satisfying Xji^v^S,) ~ 2j — 1. Its eigenfunctions 
{VjiC^ ')) j>i fo^™ orthonormal basis of L^(IR), and (pj{^, •) ~ /ij as — )• in 
L norm. As a consequence A and T have the same domain and D{Ap/'^) = W. 
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We will prove these facts for the particular class of potentials we will consider 
(see Lemmas 16.21 and 16.31 below) . 

The parameter e > will be small so that we can apply Theorem 12.31 and 
v > will be small too, so that we have a suitable perturbation theory for the 
operator A. We now explain the restriction z/ > CqE. The aim of this section 
is to construct a potential V so that Theorem 12.31 applies, and in particular, 
(|2.3p has to be satisfied. Small values of k, I in ()2.3p and the asymptotics of 
Lemma [631 give Cu > a. This together with the condition e < eoain Theorem 
12.31 vields the result. 

We fix a finite subset ^ of N of cardinal n. Without loss of generality and 
in order to simplify the presentation, we assume J7 = {!,••• ,n}. We then 
expand u and u in the basis of eigenfunctions using the phase space structure 
of the introduction, namely we write 

n 
n 

i=i i>i 

where {8, y, z, z) G VP = T'^ x W x x (recall that is the space £% with 

^(j) = j^^'^) regarded as variables and / E are regarded as parame- 
ters (here M+ denotes the set of non negative real numbers). In this setting 
equation (j6.ip reads as the Hamilton equations associated to the Hamiltonian 
function H = N + P where 

n 

N = ^Xj{iy,^)yj + ^Aj{u,OzjZj, 
i=i i>i 
Aj(i/,0 = Aj+„(i^,0, G{u,u) = (u'u)'^+i and 
(6.2) 




^iVj + Ij)''^ *^''/3i(C,x) + ^Zj(^j+„(C,x))dx. 
i=i i>i 

For the sequel we fix {Ij)i<j<n- We assume that {9, y, z, z) G -D(s, r) for some 
fixed s,r > (recall the definition ()2.4p of D{s,r)). There is no particular 
smallness assumption on s, r, we only have to take r > with r < mini<j<„ Ij 
so that {yj + Ij)^^'^ is well-defined. 
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We now show that we can construct a class of potentials V so that Theorem 
2.31 applies. 



6.1. Definition of the family of potentials V. — 

Let {fj)i<j<n be the dual basis of {h'j)i<j<n, i-e. (fj) € Span]g(/if , . . . , /i^) 
and Jjg = 6jk for all I < j, k < n. 

We say that a = {ak)k>n+i S if — ^ < Ofc < ^ for all A; > n + 1. We endow 
the set of such sequences by the probability measure defined as the infinite 
product (/c > n + 1) of the Lebesgue measure on [—1/2, 1/2]. Then define 

9{x) = ^ afce"^/i2fc-i(\/2x), 

k>n+l 

and for ^ = (6, . . . , Cn) G H = [-1, 1]" and 

n 

(6.3) V{^,x) = Y,Uk{x) + ^ig{x). 

k=l 

The spectral data ipj and Xj are defined by the spectral equation 

(6.4) {-dl + x^ + uV{C,x))^j{C,x) = X,{uOVj{C,x), 

and we assume that the (ipj) are L^— normalised (||v'j(C) = 1 for all ^ G IT 
and j >!)■ Moreover, in order to define (pj uniquely, we impose {{pj,hj) > 0. 

In the sequel we need a particular case of estimates proved by K. Yajima &; 
G. Zhang 



Lemma 6.1 ([22]). — For all 2 < p < oo there exists q > and C > so 
that for a// ^ G n and j > 1 

(6.5) ii<^i(c,-)iiLP(R) <cr". 

The next result is the key estimate in our perturbation theory. 

Lemma 6.2. — There exist q > and C > so that for all £ U, v > 
and j > 1 

(6.6) \\ipji(,-)-Vj{ri,-)\\v^<C,^\C-r]\j-''. 

In particular ||v5j(C;') ~ ^jllL2 < C'^^ICIi""; which shows that the tpj are 
close to the Hermite functions in norm. 

Proof. — In the sequel, we write (pj{^) instead of ipj{^,-). For ^, 77 E H, we 
compute 

A{u^)ipj{r]) = {-dl+x^+i^V{tx))^j{v) = >^jit'V)^jiri)+HV{tx)-Virj,x))^j 
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Thus by (j6.3p and (j6.5p there exists a > such that 

< ^\\y{o-vmLAWM\\L^ 

(6.7) < Cz.|e-r?|r". 

Choosing r/ = in (|6.7p . and as (/?j(0) = hj and Aj(0) = 2j — 1, we get 



< Cvr^\\{A{vi)-{2j-l)) 

The previous estimate together with the general formula which holds for any 
self-adjoint operator || (^(i/^) — (2j — 1)) H^^a.^^.^ = dist(2j — 1, (t(^(i/ ^))) 
gives dist(2j — \,a{A{v < Cvj"", where a{A{i'^)) denotes the spec- 
trum of A^i^S^). A similar argument, taking ^ = in (|6.7p . leads to 
dist{Xj{uri),a{T)) < Cvj-"'. Thus for all j > 1 

(6.8) Xj{uO = 2j-l + uO{r''). 

Using that {ipk{0)k>i is a Hilbertian basis of L^(R), we deduce 

Iki(^) - {Vji0,v^jiv))¥^ji0\\l2 = ^K'^jiv) - {vj{O.Vj{'n))^j{€),'^k{i))? 

k>l 

(6.9) = Yl \{vM,v^km\'- 

With the same decomposition, we can also write 

||(A(z.0-A,-(zvr?))(^,(7?)||i, = Y,\{{AiuC)-\j{i^v))Vjiv),Mm' 



k>l 



k>l 



k>l 

(6.10) > Yl \MO,Vj{ri))\^ 

k>l,kf^j 

because by ()6.8p | Afc(z^ ^) — Aj(z/ 7/)| > 1 for k ^ j uniformly in ,^,r/ and uni- 
formly in v small enough. Now by ()6.7p . ()6.9p and (|6.10p we deduce that 

W'Pjiv) - {^j{0,'Pj(.ri))^j{0\\l2 < Ciy\^-r]\j-''. 

In particular, taking the scalar product of ipj{r]) with — {£,)■, ifj{r]))ifj{^), 
we obtain 

|i-((^,(e),<^,(7?))2| <cHc-r?ir". 
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The last two estimates imply H^JjlO ~^ji''l)\\L^ < Cz^j^ — ry|j " which was the 
claim. □ 

Lemma 6.3. — We have the following asymptotics when v — \ 
(6.11) \{yi) = 2j-\-,vi,-,o{y), Vl<j<n, 

(6.12) 

n „ 

k,(yi) = Xj^nivi) = 2(i + n)-l + i/ / (/fe+<^ifc<?)/i'+, +o(i/), Vj > 1. 

fc=i 

Proof. — We first prove ()6.1ip . We differentiate equation (j6.4p in 

take the scalar product with ^j{C) and the selfadjointness of A[u ^) gives 

(6.13) ^^dir^^'h^'^^'"'^''''^^^- 



Now by ((621) 

ifk + Sikg){^]{0 -h^j)\ < \\fk + Sikgh^ yjiO + hjh^yjiO - hj\\L2 

< C7||v^,(0-/i,||l2 ^0 

when V — > 0. Thus by definition of the fk and g and by estimate (j6.13p . we 
obtain that for all 1 < j < n 

n „ 

A,- (I/O = 2j - 1 + 1/ V^fc / (/fe + 5ikg)h] + o{v) 

k=i -^^ 
= 2j -l + u^j + o{u), 

which is ()6.1ip . 

The asymptotic of ()6.12p is proved in the same way. Observe that we can 
prove a better estimate on the error term using ()6.8p . but we do not need it 
here. □ 

6.2. Verification of Assumptions [1] and [2l — 



Lemma 6.4. — There exists a null measure set M C Zn such that for all 
a G Zn\M we have for all 1 < p,q, with p ^ q 

(6.14) [ (f,+g)hl^^^Z, 
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and 

(6.15) / {h+g){hl_,p±hl^^)^Z. 

Jm. 

Proof. — For j > 1, the Hermite function hj reads hj{x) = Pj{x)e~^^^'^, where 
Pj is a polynomial of degree exactly (j — 1), and Pj is even (resp. odd) when 
{j — 1) is even (resp. odd). We have Span^{hi, . . . ,hn) = e~^^/^M„_i[X]. 
Thus we deduce that there exist (nkj) so that 

j 

(6.16) hjix) = ^^fcj/i2fc-i(\/22;), 

fc=i 

with /ijj 7^ 0. 

We assume that q < p. The application 

(a„,a„+i,...) I — > / {fi+ g){hl^p±hl^g) 

JM. 

is a linear form. In order to prove (|6.15p . it suffices to check that this linear 
form is nontrivial. According to ()6.16p and to the definition of /i and g, the 
coefficient of On+p is 

e-^^'+P^ fin+p,n+p [ /i^(„+p)_i(\/2x)dx = e-("+P)/i„+p,„+p/^/2 ^ 0. 

Therefore for fixed p, q (j6.15p is satisfied on the complementary of a null mea- 
sure set Mp,q- Finally, (j6.15p is satisfied on Zn\M where M = ^p,q>i-^p,q- The 
proof of ()6.14p is similar. □ 

In the sequel we fix a G Zr\M so that Lemma 16.41 holds true. We are now 
able to show that Assumption [1] is satisfied. Recall that in our setting, the 
internal frequencies are A(i^^) = (Aj(i^^))i<j<„ and the external frequencies 
are = {Kj{vi))j>i with kj{v^) = Xn+jH)- 

Lemma 6.5. — There exists i'q > so that for all < i' < i^q we have 

(6.17) Meas(^{C(^n : k-X{uO + l-A{uC) = 0}^ =0, ^ {k,l) e Z, 
and for a// ^ G 11 

(6.18) /•A(i/O/0, Vl<|/|<2. 

Proof. — We prove (j6.17p by contradiction. Let {k,l) £ Z. In the case \l\ = 2 
in (j6.17p we can write 

n 

k-X{uO + l- A(l/e) = Yl ^M^O + ^n+p{vi) - Xn+q{yi) := F(l/0, 
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for some p^q > 1. Now if ()6.17p does not hold, F : — > M is a function 
which vanishes on a set of positive measure in any neighbourhood of 0, thus 
F(0) = and for all 1 < /c < n, ^■(0) = 0. By Lemma 16.31 these conditions 
read 

n 

^(2j - l)kj +2{p-q) = and 

(6.19) k,+ [{f,+6,jg)ihl+p-hl^^) = 0, Vl<j<n. 

In particular for j = 1, ()6.19p is in contradiction with ()6.15p . 
The case |/| = 1 is similar, using ()6.14p . 

It remains to prove ()6.18p . For all j > 1, Aj(i/^) — > 2j — 1 when v — > 0. 
Hence ()6.18p holds true if v is small enough. □ 



We now check Assumption [21 Firstly, thanks to (|6.8p we have that for 
j,k > 1, \Aj{uC) - Afc(i/OI > \j - k\ and > j- Then by and 



< I'lC - V\\\fk + '^ifeS'll r2 sup \\ipj+n{C, 

sen 



and Assumption [2] is fulfilled. 



6.3. Verification of Assumptions [3] and (4]. — 

Recall that for p > 0, W = D{Tp/'^) is the Sobolev space based on the har- 
monic oscillator. Thanks to ()6.6p and (|6.8p . we also have W = D{Ap/^{uO) 
for all > small enough and ^ G IT. Observe that Ti^ is an algebra and the 
Sobolev embeddings which hold for the usual Sobolev space are also true 
here, since C Ti^. 

Let u = X]j>i Cij^j- Then u G Ti^ if and only if aj G £p. 

We now check the smoothness of P and the decay of the vector field Xp. 
Let p > 2 so that we are in the framework of Theorem 12.31 Since G{u, u) = 
{uu)"^~^^ in ()6.2p . we have 



(6.20) 



P = e [ 

Jr 



44 



BENOIT GREBERT & LAURENT THOMANN 



We first show that §§- ^ it We have 



(6.21) 



— = e[m + 1) / (fj+nU u 
Czj Jr 



thus ^ is (up to a constant factor) the (j + n)th coefficient of the decompo- 
sition of ■u™'u™^^, and this latter term is in Ti^ (because T-L^ is an algebra), 
hence the result. The other components of Xp can be handled in the same 
way, and we get Xp £ . 
By ()6.20p and Sobolev embeddings 



(6.22) 



I Dl ^ II ||2(m+l) . II ||2(m+l) 

sup \P\ < e\\u\\ '^i^ < e\\u\\^^ '. 

D(s,r)xn 



Similarly, using ()6.2ip and 
dP 



dP 

9yj 



ei{m + l){yj + 
^^—iVj + Ij) ' 



,m+l —m 



it is easy to see that sup^)^^ ,,)xn \^p\r < Ce. We now turn to the Lipschitz 
norms. Let ^, ry G 11 

m)-Piv)\ < Ce||n(0-n(r?)||i.(h(0||S^i+i) + ll^(^)llSti+i)) 

(6.23) < Ce\\u{0-u{r,)h2\\u\\'^,-^\ 

Now by dSSl) 



\H0-u{v)\\l2 < Cj2\\fjiO -Vj{v)\\L'^ + ^f\zj\\\fj+n{0 - fj+n[ri)\\L2 

j=i i>i 
(6.24) < C\^-vl 

where in the last line we used Cauchy-Schwarz and the fact that (-Zj)j>i G Ip 
with p > 2. Then ()6.23p and ()6.24p show the Lipschitz regularity of P. We 
can proceed similarly for Xp. 

It remains to prove the decay estimates of Assumption HI Using ()6.2ip , ()6.5p 
and the Sobolev embeddings, we obtain 

dP 

< e{m + l)\\ipj+n\\L°°(R)\\ufj2'^+i < Cer'^\\u\\:^;' % 



dzj 



and similarly, from 

d^p 



dzjdzi 



em{m + 1) / <fj+n'fi+nU 
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we deduce 



dzjdzi 



2m, ^ ^ ^j\—ce\\ „,\\2m 



< Ce\\(fj+n\\L°°{R)\\fl+n\\L°°{M.)\\u\\]%z < £C{jl) 



U 



HP- 



The estimates of the Lipschitz norms are obtained as in ()6.23p . (|6.24p and 
using ()6.6p . 

As a conclusion Assumptions [1] - 2] are satisfied and we can apply Theorem 
2.31 with some /3 > if e > is small enough. Recall that 11 = [—1, 1]". 



Theorem 6.6. — Let m > 1 and n > 1 be two integers. Let V{£,, •) be the n 
parameters family of potentials defined by (j6.3p . There exist eq > vq > 0, 
Co > and, for each e < Eq, a Cantor set Ilg C 11 o/ asymptotic full measure 
when e — )• 0, such that for each ^ G He and for each C^e <v<V{), the solution 
of 

(6.25) idtu + dlu- x^u- vV{i,x)u = e\u\^'^u, (t,2;)GMxM 
with initial datum 

n 

(6.26) uo{x) = Y,iy'e''^^,{tx), 

i=i 

with (/i, • • • , In) C (0, 1]" and 6 £ T", is quasi periodic with a quasi period uj* 
close to uq = (2j — l)j=i- |w* — Wo I < Cu. 

More precisely, when 9 covers T", the set of solutions of ()6.25p with initial 
datum ()6.26p covers a n dimensional torus which is invariant by ()6.25p . Fur- 
thermore this torus is linearly stable. 

Remark 6. 7. — From the proof it is clear that our result also applies to any 
non linearity which is a linear combination of |up"^it. Moreover, under ad hoc 
conditions on the derivatives of G, we can admit some non linearities of the 
form ^{x,u,u) (i.e. depending on x) in ()6.ip . Also we can replace the set 
{1, • • • , n} by any finite set of N of cardinality n. 



7. Application to the linear Schrodinger equation 

In this section we prove Theorem 11.21 following the scheme developed by H. 
Eliasson and S. Kuksin in [7] for the linear Schrodinger equation on the torus 
with quasi-periodic in time potentials. 

The setting differs slightly from Section [U] since now we are not considering 
a perturbation around a finite dimensional torus but we want to construct a 
linear change of variable defined on all the phase space. Consider the equation 

(7.1) idtu = —d'^u + x'^u + eV{tuj,x)u 
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where V satisfies the condition (jl.Sp . Recah the definition of the phase space 
pp = j"^ X X £p X £p. Recall also that hj, j >l denote the eigenfunctions 
of the quantum harmonic oscillator T = —5^ + and that we have Thj = 
(2j — l)hj, j > 1. Expanding u and u on the Hermite basis, u = Ylj>i ^j^j^ 
u = X]j>i equation ()7.ip reads as a non autonomous Hamiltonian system 



(7.2) 
where 



J Zj = -i{2j - l)zj - ie-^Q{t, z,z), j > 1 
\ % = «(2j - l)zj + ie-£-Q{t, z, z), j >l 

Q{t,z,z)= / V{ujt,x)[y Zjhj{x))[y Zjhj{x))dx 



anc£f3 {z, z) ^ (.\x l\. We then re-interpret (|7.2p as an autonomous Hamilto- 
nian system in an extended phase space 



(7.3) 



where 



Zj = -iC^j - 'i-)zj - ie-§-Q{e, z, z) j > 1 
% = i{2j - l)zj + ie£-Qi9, z, z) j > 1 
Qj = ujj j = 1,... ,n 



Q{9,z,z)= j v{e,x){Y,zM^)){Y.^M^)y^ 

j>l i>l 

is quadratic in {z,z). We notice that the first three equations of ()7.3p are 
independent of y and are equivalent to ()7.2p . Furthermore ()7.3p reads as the 
Hamiltonian equations associated with the Hamiltonian function H = N + Q 
where 

n 

N{cj) = Y^^jyj + Y.{2j - \)zfz,. 
i=i i>i 

Here the external parameters are directly the frequencies w = (wj)i<j<n £ 
[0, 27r)"' =: H and the normal frequencies = 2j — 1 are constant. 

7.1. Statement of the results and proof. — 

Theorem 7.1. — There exists £o > such that if < e < £q, there exist 
(i) a Cantor set H^ C H with Meas{Ii\I\.e) — )• as e — ;> ; 



'*^For the moment we work in 1% x £2, the largest phase space in which our abstract result 
applies. 
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(it) a Lip schitz family of real analytic, symplectic and linear coordinate trans- 
formation $ : He X —7- V'^ of the form 

(7.4) <^Uy, 0, z) = {y + ^z- M^{e)z, e, L^{e)z) 

where Z = {z^z), L^[9) and M^^{6) are linear hounded operators from 
ip X into itself for all p > and L^{9) is invertible ; 
(Hi) a Lipschitz family of new normal forms 

n 

such that on x 

Ho<l> = N*. 

Moreover the new external frequencies are close to the original ones 

- ^\2(5,n, < ce, 

and the new frequencies satisfy a non resonant condition, there exists a > 
such that for all G He 

\k-uj + l-^*(uj)\ > a , , (k,l)eZ. 

Notice that in the new coordinates, {y',9\z',z') = ^^^{y,6, z,z), the dy- 
namic is hnear with y' invariant : 



' i', = in*z'. j > 1 



(7.5) 



^ = -inp' j>i 



3 
:l 

'j """J^J 



As (|7.ip is equivalent to (|7.3p . this theorem imphes Theorem ll.2l In particular 
the solutions u{t,x) of (|7.ip with initial datum uo{x) = l^j>i Zj{0)hj{x) read 
u{t,x) = ^j->iZj{t)hj{x) with 

{z,z){t) = L^{Lot){z'{Qy^*\z'{Q)e-'^*') 

and (^'(0),z'(0)) = L-\Q){z{Q),m)- 
Thus 

n(t, x) = ^ Tpji^t, x)e*'^i* 
i>i 

where i^j{e,x) = E,>i[^^(^)^jHO)(^(0), ^(0))],/i,(x). 

In particular the solutions are all almost periodic in time with a non resonant 
frequencies vector {uj^VL*). Furthermore we observe that iljj{ujt,x)e^^^^ solves 
()7.ip if and only if $7* + /c • a; is an eigenvalue of ()1.10p (with eigenfunction 
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'ipj{9,x)e^^'^). This shows that the spectrum of the Floquet operator (|1.10p 
equals {Q* + k ■ uj \ k £ I^"" , j > 1} and thus Corohary 1.4 is proved. 

Remark 7.2. — Although <I> is defined on V^, the normal forms N and N* 
are well defined on only when p > 1/2. Nevertheless their flows are well 
defined and continuous from V'^ into itself (cf. (|7.5|) ). 

Proof. — Let 11 C 11 be the subset of Diophantine vector of frequencies w, i.e. 
having the property that there exists < a < 1 such that 

(7.6) \k-uj-b\>—^, A;eZ"\{0}, 6eZ 

\ky 

for some r > n + 2. It is well known that Meas(n\n) = 0. Further this 
Diophantine condition implies that 

\k-u + l-n\> a ^ , , (k,l)eZ, 

since / • G Z and if (/) < 27r|A;| then > a^q^ while if (/) > 2TT\k\ then 

|A: • w + / • > 2{l) — 2Tr\k\ > (l) > a j:^^- Thus Assumptions [T] holds true. 
Further as the normal frequencies ilj = 2j — 1 are constant, [2] is satisfied. 
We now show that Assumption [3] holds. Because of the assumptions on the 
smoothness of V, the only condition which needs some care is that {§^)k>i £ 
ip. We have 



= [ V{6,x)hkudx, 



which is the kth coefficient of the decomposition of V{9, x)u in the Hermite 
basis. Thus {§^)k>i G ^^'^ '^^^y V{6,x)u G T-L^ which is true since 
uGH^ and V and d^V are bounded. 

We turn to Assumption [H Recall that by (j6.5p . for all 2 < r < +cxd, there 
exists /3 > so that ||Lr-(iR) < Cj~^. On the other hand, by assumption V 
is real analytic in 9 and L"? in x for some 1 < q < +oo. Consider 1 < ^ < +oo 
so that i + = = 1, then with Holder, we compute 



dQ 



V{9,x)hkudx\ < sup \\V{9,-)\\Lg(js,)\\hku\\L<im 

6'G[0,27r]" 

< sup ||T^(6', •)||L9(R)l|/ife||L29(M)lkllL29(IR) 
6'e[0,27r]" 



< Ck 
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Similarly. 



dzk&zi 



V{e,x)hkhi(lx\ < sup 11^(6*, •)||L9(IR)ll^fc|lL29(M)||/ii|lL29(M) 
6'e[0,27r]" 



Therefore, Theorem 12.31 applies (with p = 2) and we almost obtain the con- 
clusions of Theorem 17.11 Indeed, comparing with Theorem \2.3\ we have to 
prove: 

(i) the symplectic coordinate transformation $ is quadratic (and thus it is 
defined on the whole phase space) and have the specific form ()7.4p ; 

(ii) the new normal form still have the same frequencies vector u ; 

(iii) the new Hamiltonian reduces to the new normal form, i.e. R* = ; 

(iv) the symplectic coordinate transformation <I>, which is defined by Theorem 
Oon each V^, extends to = T" x M" x x ij. 

Actually, at the principle Q is homogeneous of degree 2 in Z and independent 
of y and the same is true for F the solution of the first homological equation 

{F,N} + N = eQ. 

As a first consequence, does not contain linear terms in y and thus lo re- 
mains unchanged by the first iterative step (cf. (|2.1ip ). Now going to Lemma 
13.51 we notice that following notations ()3.34p . bo = bi = a = 0. Therefore 6 
remains unchanged {9 = 0) and the equation for Z reads Z = JA(9)Z which 
leads to Z{t) = e^-^^(^)Z(0) (see dOOjl ). Thus Z(l) = L^^\e)Z{0) where 

(9) = e"^^^^^ is invertible from onto itself. 
In the same way, y(r) = —^VgA{9)Z{T) ■ Z{t) (see p.46p ) which leads to 
y(l) = y(0) + iZ(0) • M^{d)Z{Q) for some linear operator M^{9). Finally the 

new error term (cf. (|2.12p ) (5+ = / {Q{t), F} o Xp dt is still homogeneous of 

^0 

degree 2 in Z and independent of y. Thus properties (i), (ii) are satisfied after 
the first step and the new error term conserves the same form. Therefore we 
can iterate the process and the limiting transformation $ = <I)^ o $^ o • • • also 
satisfies (i) and (ii). Furthermore the transformed Hamiltonian as well as the 
original one is linear in y and quadratic in Z and thus (iii) holds true. 
It remains to check (iv). This follows from the fact that <i> is a linear sym- 
plectomorphism and thus, as remarked in |12^ Proposition 1.3'], extends by 
duality on £p x £p for all p S [—2, 2] and in particular for p = 0. □ 

Proof of Corollary \1.3[ — The point is that, when V is smooth with bounded 
derivatives, the perturbation Q satisfies Assumption 3 for all p > 0. That 
is Xq maps smoothly into itself. Therefore Theorem 12.31 applies for all 
p > 2 and by ()2.6p . the canonical transformation $ is close to the identity in 
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the T^P-norm. Since in the new variables, (y' , 9' , z' , z') = <I> ^{y,9, z, z), the 
modulus of Zj is invariant, we deduce that there exist a constant C such that 

{1-Ce)\\zmp<\\z{t%<{l + Ce)\\zmp 
which in turn implies 

(1 - eC)\\uo\\nv < \\u{t)\\nP < (1 + eC)\\uo\\nP, Vt e M. 

□ 

7.2. An explicit example. — 

Consider the linear equation 

(7.7) idtu = -dlu + x^u + eV{toj)u 

where V : T" — > M is real analytic and independent of x G M. Up to a 
translation of the spectrum, we can assume that V{0) = 0. Notice that this 
case is not in the scope of Theorem 17. 11 since V does not satisfy (|1.8p . 
We suppose moreover that V = and that w G [0, 27r)" is Diophantine 
(see fm ). 

Define v{t,x) = e''^ ^o^^'^''^'^'u{t,x). The function u satisfies 1^ iff v 
satisfies idtv = —d'^v + x'^v. This latter equation is explicitly solvable using 
the Hermite basis, and the solution of ()7.7p with initial condition uo{x) = 
Yl'jLi otjhjix) then reads 



u{t, x) = e'^^o ^ ajhj{x)e 



i{2j-l)t 



Write V{e) = J2 

ake^ ' . Then, as uj is Diophantine, we can 

compute / V{us)ds = -i Yl ■^{e'^''^^ - 1), and W defined by 
W{6) = exp (e ^ ^ (e^^'^ — 1)) is a periodic and analytic function in 

oo 

0. Finally, u{t,x) = ajVF(a;t)/ij(x)e*''^-'~^''* is an almost periodic function 

in time (as an infinite sum of quasi-periodic functions). 

We can explicitly compute the transformation <I> in ()7.4p . Here the Hamil- 
tonian reads H = N + Q with Q = V{9)Yjk>i \^k? ■ Set ^{y' ,9' , z' ,z') = 
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(y, 9, z, z) where 

' zj = w{e)z'^, z,=WWjzj, j>i 

k&^,k^O l>l 

Then a straightforward computation gives 

n 

Ho^y',9',z',z') = + ^(2i - 

Therefore in this case = 2j — 1. 

Finally we study the spectrum of the Floquet operator associated to the 
equation ([TT]) . Observe that W{iot)hj{x)e'^'^^~'^^* solves ([ZZD if and only if 
any 2j — 1 + k ■ uj (with j > 1 and k £ Z"') is an eigenvalue of (jl.lOp (with 
eigenfunction W{6)hj{x)e^^'''). This shows that the Floquet spectrum is pure 
point, since linear combinations ofW{9)hj{x)e^^'^ are dense in L^(M)(X'L^(T"). 



A 

Appendix 

We show here how we can construct periodic solutions to the equation 

j idtu + dlu-x^u = \u\P~^u, p>l {t,x)€RxR, 
\u{0,x) = f{x), 

thanks to variational methods. This is classical, see e.g. .20. and for more 
details. See also [3]. Recall that for s > we have defined the Sobolev space 
T-L^{M.) = L'(T^/^), where T = —d^ + x'^ is the harmonic oscillator. We also 
define ?^°^(M) = ns>o?^''(M). We then have the following result. 

Proposition A.l. — Let fj, > 0. Then there exists an L'^(R.)- orthogonal fam- 
ily if^)j>i S ^°°(M) with llv''' ||l2(k) = A* o.'iT'd a sequence of positives numbers 
{Xj)j>i so that for all j > 1, u{t,x) = e~^'^^^ ip^ {x) is a solution of (|A.ip . 

Proof. — We look for a solution of (|A.ip of the form u{t,x) = e~^^''ip{x), 
hence ip has to satisfy 

(A.2) (_92+^2)^^^^_|^|p-l^_ 

Let /i > 0, denote by Efj_ the set 

E^ = ^ip£ n^{R), s.t. \\(p\\l^r) = /^}, 
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and define the functional 
which is on E^. 

Then the problem min J{f) admits a solution Lp^ , and Lp^ solves ()A.2p for 
some A = Ai > 0. 

Indeed, by Rellich's theorem (see e.g. page 247]), for all C > 0, the set 
jv? G s.t. ||9?||l2(r) = 

is compact in L^(M) (observe that we have used the Sobolev embedding T-L^ C 
L'P^^ which holds for any p > 1). Then, if tpn is a minimising sequence of J, 
up to a sub-sequence, we can assume that — > G in -Z>^(M). Finally, 
the lower semicontinuity of J ensures that (f^ is a minimum of J in i?^, and 
the claim follows. Moreover, Ai is given by 



Ai = 

Now we define the set Ej^ = E^^ D {{f,f^)L^[R) = O}. Similarly, we may 
construct c/?^ G Ej. so that J{(p'^) = min J((p)- The orthogonality condition 

implies in particular that t/?^ 7^ (/j-*^. Let A; > 1, and assume that we have 
constructed {^)i<j<k so that ((/?*, (/J-^)/^ 2 = /x^^jj for all 1 < i,j < k. Define 
the set 

E^^ = E^n{{^,^)L2=0, i<j<k}. 
By Rellich's theorem, the set 

jv? G n^{R), s.t. ||v3||l2(i;) = 

J liid.ipf + xV') + ^Iv^r^' < (^,V^'>L2 = 0, 1 < J < fc}, 

is compact in L^(M) and we can construct 99'^+^ G so that J((/?'^^"^) = 
min J{(p). Then ip^~^^ is a nontrivial solution of (|A.2p with 

A,+i = i / (a./+i)2 + x2(/+i)2 + 
fj- J 

The regularity c^-' G 7i°° is a direct consequence of the ellipticity of the operator 
-d'i + x^. □ 
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Remark A. 2. — Of course, the proof can be generalised to a larger class 
of nonlinearities in (jA.ip . In particular, we can deal with the nonlinearity 

— 1 P+3 

— ep|^ u with e > provided that p < 5 and that Sfi 2 > is small enough. 
Indeed in that case, thanks to the Gagliardo-Nirenberg inequality we have 



(p-l)/4 



P+l 

and the nonlinear part of the energy can be controlled by the linear part, 
which enables us the perform the same arguments as previously. 
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